Transverse dipole-dipole effective interaction for sheet arrangements 
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We have succeeded to develop a model pair interaction which when added to a system of interacting 
particles can be tuned to arrange the interacting objects into sheets. The interaction is based 
on the decomposition of the dipole-dipole interaction into two components, one parallel and one 
perpendicular to the connecting line between the dipoles, and keeping only perpendicular part here. 
Various aspects of this simple interaction are discussed, in particular in connection to two recent 
papers on self assembly of carbon nanostructures. On the other hand, the features discussed are 
quite general and might be of interest in different areas of microscopic modeling. 
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In many areas of physics there is often a need for simple 
effective interactions. Perhaps the most famous might be 
the Lennard- Jones potential or the spin-spin interaction 
of the Hcisenberg ferromagnet, two models which are well 
known and used far outside the original applications. Our 
aim has been to understand the emergence of various geo- 
metrical arrangements and we have succeeded to develop 
a model two-body interaction which when added to a 
system of interacting particles would arrange them into 
sheets instead of the expected three-dimensional struc- 
tures. This model interaction can primarily be of use 
in atomistic simulations, but possibly also for systems 
consisting of more complex particles, in principle at any 
scale. 

We start from the usual form of the dipole-dipole in- 
teraction 

U{ri2) = — 3 [mi • m2 - 3 (mi ■ ei2) (m2 • 612)] (1) 

where ri2 = ^12612 is the vector connecting the two 
dipoles and 612 its unit vector. It has probably been no- 
ticed by many, but apparently not discussed by anybody 
in printed form, that eq. [l]can be rewritten as 



U ~ /(ri2) m;^ • mj — 2m^ • mj 



(2) 



The notation used follows from decomposition of both 
dipoles into two components 



m^ei2 + 
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The figure[T]reminds us of the situations when each of the 
terms dominate. For the real dipoles the two terms are 
added exactly in the given radial form, but effective in- 
teractions may naturally assume any radial form and any 
relative strength of the two terms suitable for the physical 
model in question. Our model consists of only the part 
containing the perpendicular terms, or alignment of the 
intrinsic vectors in direction perpendicular to the con- 
necting line of the two objects. We propose to call this 
model interaction "transverse dipole-dipole interaction", 
or TDDI. The effective pair potential is chosen as 




FIG. 1: The dipoles attracting each other when positive max- 
imum scalar product (a) due to the —2 term, and when the 
scalar product is negative (b) due to the first term. There is 
a zero potential and thus zero force (c) when the dipols are 
forced to remain perpendicular. 



where all the vectors e.y, m^ 



m," are unit vectors and 



5N)G(1 



I) (3) 



the function g{rij) is an arbitrary function suitable for 
the model system in question. It can be of short range or 
long range, but preferably reaching at least two nearest 
neighbours to define the plane. The argument of func- 
tion G{u) is zero when the two vectors point in parallel 
lines, in both directions. Thus G{u) must be sufficiently 
strongly peaked at w = 0. 

Considering first just four particles with such interac- 
tion, this term will be zero only when all four particles 
are in the same plane, having thus a minimum. If there is 
also some other additional interaction leading to equilib- 
rium distances, both the discussed four particle as well as 
more particle systems would prefer sheets or single lay- 
ers, much like the graphene sheets, as long as the range 
of g{ri2) covers at least two neighbours. If there is an in- 
teraction arranging the objects at certain distances and 
perhaps also with well defined bonding angles, but with- 
out the ability to assure a planar arrangement, adding 
a Active model vector - transverse dipole (TD) - to ev- 
ery particle, the TDDI would lead to sheet-like behavior. 
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even if it is only a two-body interaction. Since paral- 
lel TDs are preferred, any three neighbours would align 
them perpendicular to the plane in which they are lying. 
When all these planes become one plane, the TDDI is 
minimized. Thus TDDI added to a suitable two or three- 
body interaction would lead to self assembly of sheets. 

Combining TDDI with additional terms is important 
because when TDDI would be the only interaction, its 
radial dependence g{r) would need to have a minimum 
to assure some equilibrium distances. This could be in 
conflict with the aim to establish a stable sheet struc- 
ture, since it could lead to a close packing situation. All 
the effective internal vectors could be partially aligned 
but a closed packed structure will result from the exis- 
tence of the minimum in the function g{r). Thus the dis- 
cussed model interaction should only be used in addition 
to other inter-particle forces. In general, the g{r) func- 
tion should be monotonically decreasing with distance, 
in most cases simply providing a cut-off. 

The proposed TDDI is a pair interaction. The limi- 
tation to the perpendicular term for a single pair leads 
only to reorientation of the two internal spins with no 
influence on the pair of particles. The same is true for 
a group of three particles. All their connecting lines are 
in one plane, so that the internal vectors can reorient 
themselves without any influence on the positions of the 
particles. When a fourth particle is added, the torque on 
the internal vectors will start forcing the particles into 
one plane. 
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FIG. 2; Eight particles, Morse interaction. Above: starting 
configuration Below: after some time, with friction 

The present paper is to some degree inspired by the 
work of Rechtsman, Stillinger and Torquato who 
show how to build structures of diamond and wurtzite 
type using special pair interactions. These isotropic pair 
interactions have minima where the radial distribution 
functions (RDF) of the target structures have peaks. The 
principle proposed here is somewhat similar. 

Tewary and Yang in ref [H construct a parametric po- 
tential for a graphene sheet. They wish to achive similar 
effect as described here, but they simulate it by combin- 
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FIG. 3: Eight particles as in fig. [1 but now with the TDDI 
forces added to the same Morse interaction. Above: starting 
configuration, the same as before; Below: after some time, 
with friction 



ing three body forces in a rather complicated way. They 
use a rather complex many body potential of Tersoff type 
with an addition of three-body terms adjusting the 
planar geometry. Our proposed intrinsic vector model 
is much simpler and more flexible, and computationally 
appears as a two-body term. We suggest that a much 
simpler model than the one of ref. [ij could be based on 
adjusting the Stillinger and Weber model 0], which is a 
two plus three body interaction for tetrahedral diamond 
correlations, to the plane geometry with the addition of 
terms of the type proposed here (see also ref. @ ) 

The effect of the TDDI interaction can be analyzed 
and explored in two dimensional case much easier than 
in three dimensions. In two dimensional case the role 
of the sheet of particles is played by a string of parti- 
cles. We have implemented the TDDI into a small pro- 
gram which solves the molecular dynamics-like problem 
in two dimensions. The interactions are Morse potential, 
which is known to lead to the triangular-hexagonal closed 
packed arrangement (isotropic forces with an equilibrium 
distance). Then we added the TDDI interaction. Results 
of two runs are shown in figures [2] and [31 The variant of 
TDDI we used here is 



U (r.y,mj,mj) 



= Un 



1 - mj 



(4) 



For the small number of particles the motion becomes 
quite disordered, therefore we have added a friction term 
to be able to see the "cooled" structures. Note that 
the starting situation is a regular, but with one parti- 
cle slightly displaced. Without any tiny displacement a 
regular structure will simply " breathe" . The difference 
between the calculations is clear and demonstrate the ef- 
fect. The physical dipole-dipole interaction where the 
first term dominates would not lead to any similar effect, 
but it would also lead to some anisotropy. 

The forces on the particles are given by the gradients 
of the potential energy, and it is instructive to inspect 
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the resulting formulae. The gradient with respect to the 
Ti coordinate of the transverse TDDI term ■ mj-, is 
obtained after some rather elementary algebra to be 



1 



ni:r = 



771 ' m:; 



(5) 



Clearly, this is the part important for the motion of the 
particles when the vectors mi, nij, are "frozen", or kept 
constant, the torque terms leading to their reorientation 
will be addressed below. The above gradient must be 
used in the evaluation of the forces on each of the parti- 
cles. For a pair potential 



W{ri,m„rj,inj) = F(ry) + Go 1 - {mj- ■ m 



(6) 



the force on i-th particle from j-th particle is obtained as 



■ {mj- ■ mf) 



(7) 



using the gradient formula given in eq. El 

Now we address the question of the torque on the vec- 
tors m.i due to their mutual interaction. While the above 
formulae are equally valid both in two and three dimen- 
sions, we limit the discussion here to two dimensions, 
since a full three-dimensional discussion becomes some- 
what lengthy. In two dimensions, there is only one angle 
specifying the orientation of each of the vectors. Denot- 
ing ai the orientation angle of i-th internal vector mi and 



9ij the direction of vector , 



cos{ai — aj)—coii{on — 9ij) cos(aj 



;) (8) 



following directly (also valid in three dimensions) from 



is spurious, since the TDs are not real quantities. When 
the optimal orienatation is reached, in normal Newtonian 
treatment the potential energy would be turned into some 
(mostly rotational) kinetic energy, leading to spurious vi- 
brations in the angles. Thus, it should be a part of the 
model to dispose off this energy, through a fictive fric- 
tion. We propose that this can be done by considering 
an over-damped motion of the model transvers dipoles, 
leading to first order equations, as described in detail in 
ref. Q. 

It might be of interest to consider a further generaliza- 
tion of the discussed model, by including both parts of 
the dipole interaction, as 



W{rij) = t(r„-)m^ • m:^ +p(rij)m| • m' 



(11) 



Here the radial functions might be such that they would 
lead to parallel layers, i.e. the p{Tij) should have a repul- 
sive part for small distances and a minimum at the de- 
sired distance between two layers. The part t(rjj) should 
be similar to the g{rij) of eq. [31 or somewhat modified 
to model even further features. 

An extensive discussion of the applications of the pro- 
posed type of model interactions is discussed in our above 
quoted paper Q . We would like to publish this shortened 
discussion of the model in this journal, since it might be 
useful in entirely different fields of research, or at least 
inspire similar applications, much in the same way as the 
surprising model of isotropic pair interactions leading to 
complex space patterns of diamond-like lattice discussed 
above (ref. Q) inspired our work. 



m^ • m^ = mj • m^ — m'l nij 



(9) 
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To evaluate the torque on m^ we need to find the gradient 
(or derivative in this case) 

d 

{mf- ■ mf) = cos(Q!i - 9ij) sin(aj - 6ij) (10) 

The torque on the m^ is zero when Ui ~ (^ij + ^y i-e. when 
it already is perpendicular to the vector r^j, connecting 
the two particles. In three dimensions the internal vec- 
tors are given by two angles, as well as the connecting 
position vector, but the principle remains the same and 
the evaluation is more complicated but straightforward. 

Possible applications of this model with fictive vari- 
ables must solve the problem of removing the energy as- 
sociated with the fictive degrees of freedom. In our ap- 
plication described elsewhere we use a model of over- 
damped angular motion of the internal vectors. The ex- 
citation energy resulting from the nonoptimal orientation 
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